Reactions at Polymer Interfaces: Transitions from 
Chemical to Diffusion-Control and Mixed Order 

Kinetics 
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ABSTRACT 



We study reactions between end-functionalized chains at a polymer-polymer interface. For 
small chemical reactivities (the typical case) the number of diblocks formed, 7Z t , obeys 2nd 
order chemically controlled kinetics, lZ t ~ t, until interfacial saturation. For high reactivities 
(e.g. radicals) a transition occurs at short times to 2nd order diffusion-controlled kinetics, 
with 7Z t ~ t/lnt for unentangled chains while t/kit and t 1 / 2 regimes occur for entangled 
chains. Long time kinetics are 1st order and controlled by diffusion of the more dilute species 
to the interface: 1Z t ~ t 1 ^ for unentangled cases, while !Z t ~ t 1 ^ and t 1 ^ regimes arise for 
entangled systems. The final 1st order regime is governed by center of gravity diffusion, 
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In many technological applications, interfaces in polymeric multicomponent materials are 
reinforced through chemical bonding |I[ . Techniques of reinforcing interfaces separating im- 
miscible polymer species involve attaching functional groups to a certain fraction of the bulk 
chains Under melt conditions, reactions produce copolymers (see fig. [I]) which serve 

as interfacial bridges, enhancing interfacial fracture toughness and yield stress after cooling 
U|. In commercial applications ("reactive processing") the two melts are simultaneously 
mechanically mixed. In addition to their direct reinforcing effect, the copolymer products 
also help to produce a finer blend morphology by reducing interfacial tension and preventing 
the coalescence of minority phase droplets during mixing ||. 

Due to the importance of such applications, recent theoretical and numerical research 
efforts have focused on understanding the reaction kinetics at polymer-polymer interfaces at 
a fundamental level []5],|^,[7||| , in the simplest case of end-functionalized chains (see fig. [I]). 
These studies have addressed only a small fraction of the available parameter space. They 
emphasized the limit of "infinitely" reactive groups, i. e. local chemical reactivities Q of order 
l/t a where t a is the monomer relaxation time. In fact, such Q values are realized only for 
exotic species such as radicals (or electronically excited groups as in model photophysical 
systems) whilst chemical reactivities employed in most reactive blending experiments are 
extremely small. A typical example is the reaction between carboxylic acid and epoxide ring 
groups for which Qt a ~ 10 -11 |§. In addition, with the exception of ref. |7|, previous studies 
have been confined to reactive chains very dilute in an unreactive melt background. 

In this letter we develop a theory for the kinetics of reactions between end-functionalized 
chains at a stationary interface separating immiscible melts A and B (see fig. [TJ). Our 
approach is quite general in terms of chemical reactivity, molecular weight, and number 
densities of reactive ends in each bulk, n A and . The convention throughout is that A 
reactants are more dilute, < n^. We will construct "phase diagrams," representing 
various classes of reaction kinetics, in the density-reactivity plane for both unentangled and 
entangled polymer melts. 

The principal conclusions of this study are as follows: (1) For the usual case of small 
reactivities, Qt a ~ 10~ 6 , the behavior is as follows. Chemically controlled 2nd order kinetics 
apply at short times, with number of reactions per unit area, lZ t , increasing as lZ t ~ t n^n^. 
If n A ~ , these kinetics persist until crowding of the interface by A-B diblocks strongly 
suppresses further reactions. However, for n A sufficiently smaller than , then before the 
interface gets saturated a crossover occurs to 1st order diffusion-controlled (DC) kinetics 
with !Z t ~ x(t) n A where x(t) is the rms monomer displacement in time t. (2) For the case 
of highly reactive groups, Qt a ~ 1, irrespective of density a crossover occurs from 2nd to 1st 
order kinetics at long times. In addition, a new 2nd order DC regime appears during which 
Tl t « x 4 ( y t)n A °n B 3 . 

Our aim is the reaction rate, j^Tlt- To begin, we observe that this requires calculation 
of the 2-body correlation function p s AB , proportional to the number of A-B reactive groups 
in contact at the interface: 4rT^t = Qha 3 p AB (t). Here h is the width of the interface and 
a the monomer size. Unfortunately, the exact dynamical equation for p s AB involves 3-body 
correlations. It is in fact impossible to write a closed exact equation for p s AB since corre- 
lation functions of all orders are coupled in an infinite hierarchy of dynamical equations. 
This complication which arises in all many-body reacting systems is typically resolved by 
approximating 3-body correlations in terms of lower order correlation functions ||10|| . How- 
ever, in ref. |TT] (where a general theory for interfacial reaction kinetics was developed) after 
postulating much less restrictive bounds on the magnitude of 3-body correlation functions 
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we obtained the following closed self-consistent equation for p s AB : 

p s AB (t) « rtfnf - X f dt' S\% p s AB (t') - Xn™ f dt 1 Sf\, p s AB (t') , X = Qha 3 , (1) 

Jo Jo 

where St ~ l/x 4 (t) is the probability density || that in the absence of reactions an A-B 
pair is in contact at the interface at t, given its members were in contact at the interface 
initially. St ~ ^/ x (t) is the probability density an A or B group, initially at the interface, 
returns to it after time t. (Eq. (|l|) is closed in terms of the degrees of freedom describing 
the reactive end-groups only. To achieve this, it has implicitly been assumed that reactions 
do not modify equilibrium internal polymer chain configurations [12| . There is considerable 



evidence from renormalization group analyses |13| that this approximation correctly captures 



all scaling dependencies, though prefactors are unreliable.) Eq. ([!]) is a closed system in 
terms of p s AB ; though spatial degrees of freedom in the bulk have been integrated out, the 
behavior of p s AB reflects bulk correlation functions. For example, the signature in eq. (0) of a 
react ion- induced density "hole" near the interface is domination of the n^n^ and S^ l > terms 
(these 2 terms then balance one another). On the other hand, a depletion hole in the bulk 
2-body correlation function corresponds to a balance of the n^n^ and terms [XT]. In the 



following, instead of directly solving eq. ([I]) for p s AB , we will rather motivate the solutions 
using physical arguments. The reader may then verify the results by direct substitution into 
eq. ©. 

Consider the reaction rate at short times. Initially, reactive groups are distributed as in 
equilibrium, p s AB w n° A n'^. Hence initial kinetics are 2nd order, lZ t ps Xtrf^n^. However, 
reactions may perturb equilibrium correlations. Consider times shorter than ti, the diffusion 
time corresponding to a distance of order the typical separation between the more dense B 
reactive groups, x(t{) = (n^ ) -1 / 3 . For £ < t\ reactions are due to the few isolated A-B pairs 
which happened to be initially within diffusive range, i.e. within x(t) of one another near 
the interface (see fig. [I]). After time t, the number of times such a pair has collided, Af co ii, is 
of order the number of encounters the A makes with the interface, (t/t a )(h/x(t)), multiplied 
by the probability it meets the B, a 3 /x 3 (t), during each of these encounters. Hence the total 
reaction probability, Qt a jV co a ~ Xt/x 4 (t) is increasing with time provided the dynamical 
exponent z characterizing monomer diffusive dynamics, x(t) ~ t 1 ^, is greater than or equal 
to 4 (the case z = 4 is marginal, but it can be shown to be similar to z > 4). Now for times 
shorter than r, where r is the longest polymer relaxation time, values of z are either 4 or 8 



TA\. This implies that the pair reaction probability is increasing with time and hence there 
exists a timescale obeying x 4 ^)/^ = X, after which all such A-B pairs will have reacted. 
Thus for t > t\ non-equilibrium correlations develop (p%b * s no longer equal to n^n^) 
and a depletion hole grows in the 2-body correlation function. TZt is then proportional to 
the number of pairs per unit area which diffusion could have brought together by time t, 
These are 2nd order diffusion controlled kinetics. 
For how long do these DC kinetics persist? For t > ti, since at least one B lies within the 
exploration volume of any A within x(t) of the interface, thus any such A must have reacted. 
For t > ti, a depletion hole thus grows in the A density field and the reaction rate becomes 
controlled by the diffusion of A species to the interface, lZ t ~ x(t)n A . These are 1st order 
DC kinetics. In summary, 

Tit ps Xtn°ln^ x(t)n^ (t* < t, < r, S conc ) . (2) 



3 



Notice that while the more dilute side controls the long time reaction rate, it is the more 
dense which determines the crossover time tj. Later, we will return to eq. (0) with explicit 
expressions for x(t). 

Certain assumptions have been made in deriving eq. It was implicit that t 2 < U < r 
since for t > r chain center of gravity Fickian diffusion takes over and x(t) ~ t 1/2 , i.e. z — I 
and A/" co ii no longer increases with time; hence our previous arguments are invalid unless both 
t 2 and t\ occur before r. Our arguments also assumed that t\ <U; were this not true, then 
an A group could have encountered many B's and yet remain unreacted. These constraints 
limit the validity of eq. (|2]) to sufficiently high Q and only. In the Q-n^ plane, these 
conditions are satisfied in the region marked S conc (for strongly reactive, concentrated) in 
figs. and |. 

Let us see first how the kinetics of eq. (fj) are modified when Q and have values such 
that t 2 > t\. Consider an A group initially within x(t) of the interface. By time t it collides 
with the interface of order (t/t a )(h/x(t)) times. Now since t 2 > ti, its reaction probability 
by ti must be small; this means each such group collides with many B's before it reacts. Now 
each encounter with the interface produces reaction with probability ~ n^a 3 Qt a . The net 
reaction probability Xn^t/xit) ~ t 1 ^ 1 ^ is thus increasing with time and becomes of order 
unity at t* m1 where x{t* m )/t* m = Xn^. It follows that for t > t^, all A groups within diffusive 
range of the interface will have reacted and a crossover occurs to 1st order DC kinetics: 

lit ~ Xtn%n% x(t)n% (t* 2 > U or t* 2 > r, W) (3) 

Notice the absence of 2nd order DC kinetics. In fact, one can show eq. @ is also valid if 
t 2 > t. This is because reaction probability for a single pair always stops increasing after r. 
In the Q-n^ plane of figs. and ^ the conditions under which eq. @ is valid correspond 
to region W, for "weakly reactive" (small Q). 

The remaining region of the Q-n^ plane is that labeled S dl1 ("strongly reactive, dilute") 
in figs. @ and ||, corresponding to t* 2 < r < t[. In this region all arguments leading to eq. 
(0) are still valid for t < r. But for t > r, z = 2 and reaction is no longer guaranteed when 
2 groups are within diffusive range. Thus the hole in the 2-body correlation function stops 
growing and the 2nd order DC kinetics cease. However, since t\ < r, reaction is certain when 
two coils overlap at the interface. Since by r there is no density depletion at the interface, 
thus for t > r, the reaction rate is proportional to the equilibrium number of overlapping 
coils per unit area at the interface: TZ t ~ {R^/^tn^n^, where R is the polymer coil radius. 
Repeating the argument leading to the calculation of t* m) but replacing a and h by R and 
Q by 1/r, leads to the conclusion that the reaction rate becomes limited by delivery of A 
species to the interface at a timescale t^ nR = r(i? 3 n^)~ 2 . In summary: 

K t ss Xtn%n% ^ x 4 (t)n%n% (R 4 /r)tn^n^ x(t)n% (t* 2 < r < t h S dil ) 

(4) 

This completes the derivation of our general results, eqs. (fj), @ and (f|). Let us now apply 
these to unentangled and entangled melts, respectively. 

1. Unentangled chains. For shorter chain lengths N, Rouse dynamics apply |14j], x(t) ~ 
t 1//4 for t < r fa N 2 t a . Ordinary Fickian diffusion, x(t) ~ t 1 ^ 2 , onsets for t > r. Following 
the steps described above, the phase diagram of fig. is constructed. Explicit formulae for 
timescales are obtained after inserting the appropriate x(t) forms. The only complication is 
that the short time z = 4 regime is marginal (J\f co a increases only logarithmically with time). 
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This leads to logarithmic corrections during the 2nd order x 4 (t) DC kinetics: 

TZ t « aH/[t a ]n{t/t h ))n%n% (2nd order DC) , (5) 

where th is the diffusion time corresponding to h. Similarly, the (i? 4 /r)t n^n^ regime 
is modified to TZt ~ -R 4 £/[ rm ( r /^)] n A . These logarithmic corrections also modify the 
timescales t^ti and t* mR . From the results presented here, they are easily determined by 
demanding continuity of reaction rates. 

As an example, consider a point in fig. belonging to region S conc . The kinetic sequence 

is: 

TZt ~ tn%n% (t/hit)n%n% t l/i n°£ t 1/2 r^ . (6) 



Entangled chains. For iV above the threshold iV e , entanglements onset. According to 
the "tube" or "reptation" model |TJ]] Rouse dynamics apply for t < t e = N 2 t a until x(t) 
reaches the tube diameter r P = N. 



72 



a. 



The "breathing modes" t 1 ' 8 regime follows, lasting 
until t h = N 2 t a , while x(t) ~ t 1 / 4 for t b < £ < r where r = (N 3 /N e )t a . Diffusion is Fickian 
for t > t. 

Similarly to the unentangled case one constructs the phase diagram of fig. |3|. The 3 basic 
regions S conc , S dl1 and W now develop fine structure since each characteristic timescale may 
belong to different x(t) regimes. (For example, for point X of fig. |3| one has t e < t\ < %, 
and % < ti < r.) Again, logarithmic corrections to the 2nd order DC kinetics arise during 
the two short time t 1//4 regimes, with TZt ~ a 4 t/[t a In^t/th^n^n^ during the first, and 
TZ t ~ r^t/\t\>\n(t/t\ ) )}n™n < B during the second (where r b = N^N^a). 

As an example consider point X marked in fig. |3], belonging to the S conc regime. Then 
eq. (0) implies: 

TZ t - tn%n% t 1/2 n^n™ -S {t/\nt)n^n™ -±> t 1/4 n% t 1/2 n% . (7) 



Finally, let us consider effects of interfacial saturation by products. When TZt reaches 
the value TZ CT i t = iV _1 / 2 a~ 2 , the diblocks start to feel one another laterally and therefore 
begin to stretch PJl5|. Now during the very late stages, when TZ t 3> TZ CT i t , incoming chains 
must overcome a very large free energetic barrier in order to react with the surface. It 
was shown in ref. § that this leads to a reaction rate decreasing exponentially with surface 
coverage, TZ t ~ e~ 9 ^*^ crit ^ , which in turn implies that TZ t increases logarithmically slowly 
in time, as demonstrated in ref. 0. (For coverages greater than but of order TZ CT i t one expects 
TZ t ~ f(TZ t /TZ crit ) where f(u > 1) = e~ 9u . The detailed form of the cross-over function / 
describes this intermediate regime. / is in principle available from eq. (|I|).) 

Thus, for very large N, in effect reactions cease at coverages large compared to TZ CT i t (here 
we neglect destabilization effects which may increase interfacial area). Equating TZ t to TZ cr \ t , 
we thus find that for all points well above a certain line (labeled t* m = t sat ) in figs. || and [| 
saturation occurs during the final 1st order DC kinetics regime at a time t sa t = r/ (Nn^a 3 ) 2 . 
Below that line, t S at = 1/ ' (Qn^n^ha^R) is smaller than t* m and saturation occurs before the 
crossover to the 1st order DC regime. 

We conclude with a brief discussion of our main findings. Since ordinary bimolecular rate 
constants correspond to Qt a ~ 10~ 6 , experimental systems typically lie below the t sat = t* m 
line of figs. ^| and |||. Then, apparently consistent with experiment fl6|], linear 2nd order 
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kinetics, lZ t ~ t, apply until surface crowding by products at £ sat . As a specific example, 
if all chains were functionalized, N = 300 and Qt a fh 10 -6 , then t sat would be of order 
seconds. Typically groups are even less reactive and t sat much longer. Note that if <C rig 
then even for small Qt a a transition occurs to 1st order DC kinetics. However even in 
this case these kinetics may not be observable since t* m is likely to exceed reactive blending 
experimental timescales (of order minutes). We conclude that typical reactive blending 
experiments involving end-functionalized polymers are not in the DC regime. This would 



seem to be consistent with recent experiment [[17 



Of considerable potential interest are experiments involving highly reactive radical species 
or strongly physically associating polymers, for which Qt a m 1. For these, a range of 2nd 
and 1st order DC regimes are predicted. Measuring these various kinetic regimes with 
laser-induced radicals thus offers an interesting probe of fundamentals of interfacial polymer 
dynamics. 
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Figures 



FIG. 1. End-functionalized chains reacting at an interface separating immiscible melts A, 
B containing functional groups, of local reactivity Q, whose initial densities are n°^,n^. 
For short times, reactions are confined to those groups whose exploration volumes of size 
x(t) overlap at the interface. The number of such pairs per unit area is x 4 (t)n°^n'§. 

FIG. 2. Unentangled melts, reaction kinetics "phase diagram" in Q — plane. Axes are 
logarithmic and units chosen such that t a — a — 1. Maximum possible density is n^a 3 = 
1/N (all chains functionalized). Kinetics are qualitatively distinct in the 3 regions S conc , 
S dll and W (defined by solid lines). The S conc /W boundary defines a critical density- 
dependent reactivity, Q = Q*. Region W is divided into two sub- regions by the t* m = r 
line, above (below) which t* m < r {t* m > r). The line t sat = t* m is shown for the special 
case of constant n^/n^ (the convention is < ^s)- 

FIG. 3. As fig. but for entangled melts. Regions S conc , S dll and W now develop fine struc- 
ture. In a given sub-region, relevant timescales occur within a given reptation diffusion 
regime thus defining a unique sequence of kinetic regimes. 
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Fig. 1 
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Fig. 2 
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Fig. 3 



